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ABSTRACT 



Severed topological reduction and decomposition techniques are developed to 
decrease the complexity of computing R u k{D), the source-to-K-terminal reliability 
of an acyclic directed network D with independent component failures. R s jc{D) is 
computed in 0( J F"J | jET | ) time when D can be completely reduced. When rot 
completely reduced. a graph D' remains such that R 3 ^{D) = 

M ' RsX^D i)~ R 3 k z {Rz) • • • Rsz {Dm) where M is a knov/n constant and the D-. are one 
or more separable components of D'. A simple scheme, exponential in ' — /Q ; Is 

given for computing R sK {D). When | V t - K\ | become too large, a truncated version 
of this scheme usually gives an excellent lower bound on R-jqiDi) and thus an 
excellent lower bound on R^{D). 

A program using these techniques has been coded in FORTRAN and tested on 
"complete" acyclic graphs and "street" networks with up to 100 vertices. Running 
on an IBM 3033AP under FORTRAN H (Extended), total CPU time for computing 
exact reliability is less than 3.5 seconds when j V'-A'j^lO and Fj<100. 
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I. INTRODUCTION 



Often, in the design and production of complex system such as computer 
network, good system reliability is a critical requirement. Network reliability is 
rarely treated explicitly and quantitatively, however, largely because of problems 
of computational complexity and because few efficient computer programs exist 
for computing network reliability. Only in the last few years have research and 
development in this area provided usable results. 

We wish, as efficiently as possible, to determine the network reliability for 
specific class of networks. The purpose of this thesis is to develop new reduction 
and decomposition techniques for computing network reliability and to show how 
computational complexity can be decreased by using these techniques. 

The network model which is used in this thesis may be thought of as a 
communication network with directed edges, allowing only one way communication 
from tail vertex to head vertex. Additionally, we impose the restriction that the 
network have no directed cycles. Such a model may be appropriate in hierarchical 
communication networks such as command networks. Also, such networks may 
arise as subproblems while computing the reliability of general directed networks, 
Agrawal [1962]. 

Networks can be divided into various classes, e.g., directed, undirected, cyclic, 
acyclic, etc. A large amount of literature exists on each class. (See Agrawal 
[1982], Belli and Provan [1961], Satyanarayana [1982], Hagstrom [i960] and 
Satyanarayana and Prabahkar [1978] for directed networks; see Ball and Provan 
[1981], Johnson [1982], Satyanarayana [1982], Satyanarayana and Chang [1981], 
Satyanaryana and Wood [1982], Valiant [1978], and Wood [1982] for undirected 
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networks; see Agrawal [1982] and Satyanarayana and Prabhakar [1978] for cyclic 
versus acyclic directed networks.) The specific model we will consider is the 
directed acyclic network, denoted D{V,E ) where V is the set of vertices and E is 
the set of edges. Each edge is represented by an ordered pair of vertices e ={u,v). 
The ordering implies that the flow of co mm unication can move along e from u to v 
only. 

Each edge e* in the network functions with probability Pi and fails with 
probability = l-p<. All edge failures are assumed to be independent events. For 
simplicity in the initial discussion, vertices are assumed to be completely reliable. 
However, the extension to unreliable vertices will be handled briefly. 

Vertex u can send communication to vertex v in D if there is a directed path 
of working edges from u to v. One of the vertices s is designated as the "source" 
or the "root". With respect to the source, there are three different measures of 
reliability which Eire usually studied in directed networks: 

(1) The source-to-terminal (ST) reliability. R st {D), is the probability that s can 

send communication to a specified vertex t in D. 

(2) The source-to-all-terminal (SAT) reliability, R s y(D), is the probability that s can 

send communication to all vertices of D. 

(3) The source-to-K-terminal (SKT) reliability, Rsk{D), is the probability that s can 

send communication to every vertex in a specified set K (these vertices are 
referred to as "K-vertices”) with, by convention, s zK. 

Of course, (l) and (2) are special cases of (3). We consider the general case (3) in 
this thesis. 

Generally, the of problem computing network reliability is NP-hard. Certain 
reduction and decomposition schemes exist which reduce the size of the problem 
while preserving network reliability, e.g., Agrawal[1982], Shogan[l973] and Wood 
[1982], and, for special classes of networks, it is thereby possible to compute 
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reliability in polynomial time In general, however, the solution time grows 
exponentially in the size of the network. In this thesis, we develop a method for 
computing reliability which, in general, is exponential in the size of D, but which is 
polynomial when the number of non-K-vertices, i.e., .| V—K |, is fixed. 

The rest of this thesis is outlined as follows. 

Chapter 2 contains a necessary definitions and graph theoretic concepts. 

Chapter 3 develops reduction and decomposition techniques to reduce the size 
of a network in order to decrease computation time. Some of these techniques are 
new. 

Chapter 4 shows how to compute network reliability exactly when the number 
of non-K-vertices is moderate. A truncated version of this technique is used to 
obtaun a lower bound on reliability when the number of non-K-vertices is too large. 
We demonstrate that this lower bound will usually be "good." 

Chapter 5 describes an algorithm to compute R s k{D) or find a lower bound on 
R s fc(D), and its implementation. The program is coded in FORTRAN with one main 
routine and 18 subroutines totaling about 1600 lines of code. 

Chapter 6 gives computational experience for the algorithm run on an IBM 
3033AP system /370 with the FORTRAN H (Extended) compiler. 

Chapter 7 is a conclusion and gives suggestions for further research 
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II. DEFINITIONS AND NETWORK REIJABH JTY THEORY 



In this chapter, we give the definitions and notation used throughout this 
thesis, and introduce basic network reliability theory. 

A. DEFINITIONS AND NOTATION 

A directed network D=(V,E) comprises two finite sets: V is the set of vertices 
and E is the set of edges. Each edge eeE corresponds to an ordered pair of 
vertices in the directed network. That is, e=(u,v ) where u.veV. The vertex u is 
called the tail vertex of edge e and vertex v is called the head vertex of edge e . 
For any two vertices u and v in D, if there exists an edge e =(u,v), then u and v 
are said to be adjacent and edge e is incident into v and incident out of of u . The 
indegree of a vertex v, denoted indegfa), is the number of edges incident into v 
and the outdegree of a vertex v, denoted outdeg(v), is the number of edges 
incident out of v . Two vertices are connected (or "communicate") if there exists a 

sequence of vertices and edges of the form u 0 (uq/Vi), V\. (vj.-ug) (•u m - l ,v m ) v m- 

This sequence is a path of length m. If the v 0 =v m , the path is a cycle. 

Consider a undirected graph G(V.E') formed by ignoring the ordering of the 
edges in a directed graph D{V,E). G is called the underlying graph of D. D is 
connected if its underlying graph G is connected. G is connected if there is an 
undirected path between all pairs of vertices. The vertex connectivity of D is the 
minimum number of vertices which must be deleted from G (along with adjacent 
edges, of course) in order to disconnect G or to create a trivial graph with only one 
vertex. D is separable if its vertex connectivity is one. Otherwise D is 
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nonseparable. In a separable graph, a vertex whose deletion disconnects the graph 
is called a cut-vertex. 

Merging a vertex v into another vertex w means removing all edges between v 
and w and forming a single vertex v' by coalescing v and w such that any edge 
incident into (out of) v or w is now incident into (out of) v'. Merging v and w is 
sometimes also referred to as "identifying” v with w. 

The SKT reliability of a directed network D, denoted R 3 x(D), is the probability 
that all K-vertices in D are connected by path of working edges and vertices from 
source vertex s. The probability that edge e* works is Pi = 1 -g*. For now, we 
assume that all vertices work with probability one. 

B. NETWORK THEORY ANT) SURVEY 

1. NP-rnmplptpnp<;«! 

We briefly discuss the complexity of the SKT reliability problem in 
terms of the theory of NP-completeness and in terms of practical computation. 
More detail can be found in Garey and Johnson [1979], Ball and Provan[l98l]. 

Algorithms to solve various problems can be broadly classified into two 
categories, polynomial-time and exponential-time algorithms. An algorithm is a 
polynomial-time algorithm if for a problem of size n, its running time is bounded 
by a polynomial in n. Any algorithm that is not a polynomial-time algorithm is an 
exponential-time algorithm. In combinatorics, a decision problem P is said to 
belong to the class NP if given a tentative solution, the validity of the solution can 
be checked in polynomial-time. P is "NP-complete" if it is equivalent to the 
satisfiability problem of a conjunctive normal form boolean expression. A problem 
is NP-hard if it is at least as hard as an NP-complete problem. 
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A problem is #P-complete (number-P complete) if its solution is equivalent to 
counting (not enumerating) the number of solutions to an NP-complete problem. 
Any #P-complete problem is NP-hard. 

2. Prohlpm Pnmplpvity 

Bedl and Provan [1981] show that the ST reliability problem in an acyclic 
directed network is a #P-complete problem. This implies that the SKT problem is 
also #P-complete since problem instances restricted to \K\=Z are #P-complete. 
Even finding an e-approximation to SKT reliability is ein #P-complete problem. 

In practice, these complexity results imply that a general problem of 
large size cannot be solved analytically in a reasonable amount of time. In fact, 
the time taken grows exponentially with size of the network. As will be seen, 
however, practical problems of fairly large size can be handled with the techniques 
developed here. Furthermore, even larger problems can be handled if one is 
willing to accept a "good" lower bound on exact reliability. 

3. RpH notion 

Many reduction schemes have been developed to help solve various 
network reliability problems. These reductions reduce the size of the network in 
polynomial time while preserving reliability, i.e., they create a new, smaller, 
network G' from the original network G ( G may directed or undirected and any 
number of definitions of reliability may be used.) such that R(G)~M R{G') where M 
is a known constant factor. In this way, a polynomial amount ci work is expended 
in order to save, it is hoped, an exponential amount of work during later, more 
general computations. In some special, but often practical cases, it is possible to 
completely reduce a network and effectively compute its reliability in polynomial 
time. 

The reductions describe below give a flavor for the different types of 
reductions which have been developed. Parallel and series reductions are well- 
known schemes so they will not be discussed. The reductions described assume 
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that vertices are completely reliable. However, most can be extended to handle 
unreliable vertices. 

(1) Polygon-To-Chain Reductions: A chain C in an undirected graph G is defined as 
a alternating sequence of vertices and edges beginning and ending ’.nth vertices 
such that all interned vertices have degree 2. Two parallel chains form a polygon 
and this polygon can be replaced by a chain in a reduction that preserves 
reliability in the K-terminal reliability problem, the analog to the SKT reliability 
problem in an undirected network. This reduction and its implications are 
developed by Wood [1982]. 

(2) 2-Neighbor Vertex Reduction: For the SKT reliability problem, any vertex u 
that has exactly t T wo neighbors v and w can be eliminated from the network if 
uf? K or [u.v.wloK and u /s . 

(3) Paralle 1-Chain Reduction: Suppose E^ is the set of all edges in a directed 
network D whose endpoints are v* and Vj. In D, a chain C is an alternating 
sequence of vertices and non-empty sets of edges where the internal vertices are 
2-neighbor vertices in D. Two chains C x and C z with identical endpoints are 
parallel chains. Parallel chains can be reduced to a single chain in the SKT 
reliability problem. Reductions (2) and (3) are developed by Agrawal [1932]. 

4 . De.compn.siti an 

Sometimes it is possible to divide a given graph into a number of 
subgraphs, such that each subgraph can be analyzed separately and the results 
combined to obtain the reliability of the whole graph. We briefly describe some of 
these methods based on vertex connectivity. 

(1) 2-Connected Graph Decomposition: Let G be a 2-connected graph with a pair of 
vertices [u,v] whose deletion disconnects the graph. These vertices are called a 
"separating pair." This separating pair partitions G into two subgraphs Gi and G? 
such that Gi(^G z =G, Gif)G 2 =\u,v\, Ei^<t> and E 2 It has been shown by Wood 
[1982] that, in the K-terminal reliability problem, each subgraph defines a chain 
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through a set of reliability computations on the subgraphs. The subgraphs may be 
replaced by the chains so as to preserve reliability. Agrawal [1982] has proven the 
analogous decomposition scheme for the SKT reliability problem on directed 
networks. (2) Strong-connected Decomposition: A strong component of a directed 
graph D is a maximal subgraph of D such that there exists a directed path from 
every vertex to every ether vertex in the subgraph. When computing R s k{D), each 
strong component can be considered separately. Computing each component 
RM) and multipling Rsic(Di) is whole graph reliability R s k{0). This 
decomposition has been developed by Shogan [1973] and Agrawal [1982]. 

Other decomposition schemes exist. Decomposition in 3-connected 
graphs is often possible, Rosenthal [1974]. Decomposition in 1-connected graphs 
will be discussed in chapter 3. 
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III. REDUCTION AND DECPMPQSmON 



Any algorithm to solve a general network reliability problem has a running 
time which is exponential in the size of problem. If by some method, the size of the 
problem can be reduced, immense computational savings can be made. Therefore, 
in this chapter, we introduce reduction and decomposition schemes which will 
often reduce the size of the problem while always preserving reliability. Some of 
these techniques are specific to the SKT reliability problem in acyclic directed 
networks while others sire of more general applicability. 

Reduction schemes reduce a parameter of the size of the problem, usually 
edges or vertices, in order to decrease computational complexity while preserving 
network reliability. On the other hand, decomposition schemes split the whole 
graph into subgraphs, computes subgraphs reliabilities separately and then 
combines these reliabilities in some way to compute overall network reliability. 

A. REDUCTIONS 

In this section, we discuss various reduction schemes which reduce the size of 
D and, consequently, reduce the complexity of computing R 3 k{D). Under any 
reduction, reliability of the network remains invariant up to a known, constant 
factor, i.e., R 3 k{D ) = M'R 3 lc {D') where M is the known factor, FC is the new' set of 
K-vertices, and D' is the reduced network. In the following illustrations, all vertices 
are represented as circles with K-vertices being darkened or shaded. 

The parallel and series reductions are standard, i.e., well-known and widely 
used. The three reduction schemes, 
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"indegree-one K-vertex contraction'', "non-K-vertex deletion” and "neck vertex to 
K-vertex reduction" are newly developed reductions. 

1. Parallel Rrigps RpdnrHnn 

In a directed network, two parallel edges e a =(u,v) and s b -{u,v) may 
be replaced a single edge e c =(u,v) while preserving reliability if p c = p a + Pb~Pa?b 
and M- 1. 

2. Series-Edges Reduction 

In a directed network, a pair of edges e a =(u,u) and e b =(v,w) such that 
indeg(u)=outdeg(i/) = l are series edges. If K, edges e a emd e 0 can be replaced 
by single edge e c =(u,w) while preserving reliability if p c = p a p b and M- 1. 

3. Tndpgrpp-Onp K-vprt-pv Contrari-inn 

Let Si=(u,v) be an edge in D such that veK and indeg(u)=l. In the 5KT 
reliability problem, communication from s to vertex vzK can only only occur via 
edge so this edge must work if the network is to work. Consequently, edge e* 
can be contracted such that R& ( D ) = M-R 3k ( D ') and M =p* 



where D’ = \ V — u -v r w - ei j 
w = u iju 



IC = 



K—v +ui, if K 
K-v—u+w, if i leK 



The validity of this reduction is easily shown. Let e^u.v) be the edge 
described above in the acyclic network D. Let denote the event that e* is 
working, emd let Fi denote the complementary event. Since R 3 x(D) is just a 
probability, the rules of conditional probability can be applied to obtain 

' R^D) =PiR sf c(D\F i )+q i R sK (D\F i ) 

- since R sK (D\F i )-Q 
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= PiRsAD') 

since no new paths are created by contraction. 

4. Nlnn-K-VprtpY nplptinn 

In the SKT reliability problem, an outdegree zero non-K-vertex can be 
deleted from the network without changing network reliability since such a vertex 
is irrelevant. In this case R,x(D) = M R s x(D-v ) and M- 1. 

5. Vprt-py tn K-yprt-pY RprhirHnn 

Suppose u and v are two vertices in D such that ueK, K and all 
paths from s to u include v. v is called a neck vertex. It may or may not be a 
cut-vertex. Since s must communicate with v before it communicates with u, 
R s k{D) = M R s (x+v)(D) where M- 1. Thus, the neck vertex to K-vertex reduction 
simply changes K to K+v. See Figures 3.1a and 3.ib. (All figures are given in 
Appendix A.) 

The neck vertex to K-vertex is termed a "reduction" because it reduces 
the parameter | V—K | in D. As will be seen in chapter 4. the general algorithm for 
computing R s x(D) runs in time proportional to 2^ V ~ K ^ so this is a reduction 
according to our definition. Furthermore, the neck vertex to K-vertex reduction 
will often allow other reductions or decomposition to take place. 

6. Qthpr Rprinptinns 

Numerous other reductions exist as exemplified in chapter 2, section 
B3. However, none of these reductions has been implemented because of their 
complexity, and therefore, no other reductions will be described here. 

B. DECOMPOSITION 

Sometimes it is possible to decompose a given network D into a number of 
subnetworks. Each subnetwork can then be analyzed separately and the results 
combined in some way to obtain overall network reliability. 
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Decomposition can result in substantial savings because, in general, the time 
to analyze the whole system is exponential in size of the system. Effort needed to 
compute R,k{D) using decomposition is the sum of efforts needed for each 
subproblem plus the effort, normally polynomial in the size of D, in performing the 
decomposition and recombining the result. Without decomposition, the necessary 
work is roughly the product of the effort needed to solve the individual problems. 

1. Cut-Vertex Decnmpnsitinn 

A connected, undirected, network G{V,E) has a cut-vertex v if there 
exist components (subgraphs) C x and C 2 such that C\C\C 2 =v and Ci\jC 2 = G . Any 
communication from one component to the other must pass through v. In this 
case we also say that v separates C x and C 2 . A network which has a cut-vertex is 
"separable" and one which has none is called "nonseparable." A directed network D 
is separable if its underlying graph G is separable. 

Let D have nonseparable components C l ,C 2 C Any component Q 

which has only one cut-vertex is called a pendant component. If a pendant 
component, Q, contains no K-vertices except possibly its cut-vertex, then it is 
completely irrelevant with respect to reliability since none of its edges are in any 
path from the source to a K-vertex. Thus, this component, excluding its cut-vertex, 
can be deleted from the network while preserving network reliability in the SKT 
reliability problem. This type of deletion may be recursively applied until only 
relevant components remain. 

Cut -vertices are not efficiently found in a directed network if each edge 
is represented in one direction, e.g., if all edges incident into each vertex are 
stored contiguously. However, if edges are stored both by tail vertex and head 
vertex, then this is exactly the representation of the underlying undirected 
network. Thus, an efficient algorithm fcr finding cut-vertices in an undirected 
graph may be used to find cut-vertices in D. A slight modification to the 
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DFS (depth-first search) algorithm due to Tarjan [1972] cam find all cut-vertices in 
0( |^ |) time. 

Once all cut-vertices have been found and ell irrelevant components 
removed, any cut-vertices remai ning which are not K-vertices are neck vertices 
and may be added to the set K. It is fairly obvious that each component must work 
by itseli if D is to work. Since edge reliabilities are independent it follows that 
R*x{D) = RzxiD x ) Rtf (D z ) ••• R 3 AD m ) where the are the nonseparable, 
relevant components of D. This is identically equal to a network which has the 
same components as D but whose only cut-vertex is the source. A more rigorous 
demonstration of this identity can be made by applying the decomposition scheme 
described next, "moving edges to source,” to all edges incident out of cut-vertices. 
Figure 3.2 demonstrates the results of cut-vertex decomposition applied to a 
network D with five components, two of which are irrelevant. In this case, 
Rsx(D) = R s k(D i)'R,x iD z ) 'R sK (D 5 ). 

2. Maying. lid gps In Source 

In this section, we discuss a technique which does not fit neatly under 
the rubric of decomposition or reduction since it only moves edges from one place 
to another in the network. However, by moving certain edge tails to the source s, 
new parallel reductions, series reductions or contractions may be made possible 
along with cut-vertex decompositions. Consider an edge e=(u,v) such that ueX. 
Edge e can be made incident out of any other K-vertex u'zK without affecting 
reliability as long as no cycles are created. A simple conditioning argument can be 
used to show this, Agrawal [1932]. The source is the obvious candidate to receive 
the edge tails since it has no edges incident into it and therefore no cycles cam 
possibly be created. Figure 3.3 shows an example where this operation creates a 
cut-vertex and allows a parallel reduction. 
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IV. COMPUTING SXT RELIABILITY IN AN ACYCLIC DIRECTED NETWORK 



This chapter gives a general method for computing SKT reliability in an acyclic 
directed network. The method is based on repeated computation of SAT reliability 
of certain subgraphs of D and, consequently, we first discuss SAT reliability. After 
deriving a method for computing SKT reliability, a simple lower bound becomes 
apparent and several heuristic arguments are given indicating that the bound's 
accuracy should be good. Actual computational examples are reserved for the 
next chapter. 



A. SOURCE-TO-ALL TERMINAL RELIABILITY 

Let D(V,E ) be an acyclic directed network and let qx-l~Vi be the failure 
probability for any edge s^E and let Ej be the set of edges incident into Vj. Also, 
let the vertices be numbered such that $ v l ,v z , • • • ,v n } is an an acyclic ordering of 
V. Then, SAT reliability may be expressed as 



E sV (D) - s communicates with all ve. 7— v n 

and at least one edge into v n works ) 
since there no edges incident out of v n by the acyclic ordering 



( 4 . 1 ) 



— V-v^){D ^Vn) 



- n m 

9 j*&i 



by independence 
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n 
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Another way to see this is to use the fundamental topological property that 
R a y{D) is unchanged if any edge {u,v) with uzK is replaced by an edge ( w ,v ) with 
the same reliability and i v<zK has a lower number than u in an acyclic ordering of 
D. Consequently, in the SAT problem, we may move ail edges so that they come 
out of the source just like the "moving edge to source" operation in section 3.2. 
Equation (4.1) is then trivially true for a network where all edges come out of the 
source. (See Figure 4.1). The complexity of computing R a y{D) is easily seen to be 
0( \E |) if all edges are with common head vertices are stored contiguously. 

B. SOURCE-TO-K-TERMINAL RELIABILITY 

Now we consider the SKT reliability problem in D. Let D{V,E) be a acyclic 
directed network with KqV. By simple state-space partitioning, 

RsJe{D) = P{ s communicates •with allvsk) 

+ 2 jP ( s communicates with all veF— u* but not with Vi) 

K 

+ P( s communicates with all vzV-Vi—Vj but not with Vi or vj) 

Vj ,Vj& K 

+ . . . 

+ P{ s communicates with all vzK but not with any veV—K) 

But, 

P{ s communicates with all vzV) = R s y{D ) 

P{ s communicates ■with all veV — but not within ) = R a (y- Vi )(D). P{ No edge into 
Vi works ) 

P{ s communicates with all veF— v i —Vj)=R a ^ Vi -. Jj )(D). P{ No edge into Vi or Vj 
works except an edge e =(v i ,i; J ) or e =(vj,Vi)) 
etc., 

In general, then. 
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where Ei is the set of edges going into v it Ey is the set of all edges going into or 
Vj except edges going between Vi and Vj, E iik is the set of all edges going into v^. Vj 
or v k except edges going between pairs of vertices in etc. 

The complexity of computing E^iD) is, via the above formula, exponential in 
the number of non-X-vertices. The toted number of additive terms in equation (4.2) 

n i 

is 2 Ci V ~ K ^ - where CP is the number of ways to choose i elements out of n 

i=0 

elements. Computing SKT reliability thus requires 0(|iT|2 |V, “ /i:i ) time using 
equation (4.2). Of course, if the number of non-K-vertices is fixed as |iT| changes, 
computation is effectively 0( | jET | ). 



C. A LOWER BOUND ON SKT RELIABILITY 

Since computation increases exponentially in the worst case, we may be 
satisfied with an approximation to SKT reliability which requires less time to 
compute than exact reliability. Since all terms in equation (4.2) are positive, any 
subset of these terms provides an approximation to SKT reliability which is a lower 
bound. Also, the product formed by lower bounds on independent components of a 
network yields a lower bound on overall network reliability. The subset of terms we 
choose to form this lower, bound consists of terms enumerated in the order given in 
equation (4.2) until some specified computational limi t is reached. The motivation 
for this scheme being fairly accurate is given below. 
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Suppose p x = p for ail e i e£’ and p is close to 1. Consider the first few terms of 
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Equation (4.2). 

R, V (D) = nji - q' E A ~ 1.0 

(.A VU (l ~ ?IE ‘ IHE ' IM 

Thus, the contribution of the terms to the total reliability sum diminishes 
rapidly unless there are a significant number of edges connecting non-K-vertices. 



Another reason that the lower bound should often be good is that terms with 
many vertices deleted may be identically zero. For example, consider an ST 
reliability problem where the shortest path from s to t has L edges in it. This path 
must contain i-1 non-K-vertices and so R,(v-v i - Vj -...~v k ){D- / Vi~ v, = 0 when 

| V-[v it vj s,£ } | <1 — 1. Of course, other terms may be identically zero, too. 

For example, in any K-terminal problem, the term corresponding to D- / u i ~Vj will 
be identically zero if all paths from s tc some K-vertexi/ must include either i/, or 

Vj. 

Computational experience in Chapter V lends credence to the accuracy of this 
lower bound. 
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V. ALGORITHM AND DATA STRUCTURE 



The algorithm presented here for reliability evaluation of directed acyclic 
networks is based on the theory developed in chapter 3 and chapter 4. The 
program is intended for general use so it was coded in the widely available 
language, FORTRAN. This language is available on most computer systems and is 
still one of the most popular languages among operations research analysts. 

A. ALGORITHM 

The objective in this section is to develop an efficient algorithm for computing 
SKT reliability in a directed acyclic network. Initially, polynomial-time reduptions 
and decomposition are used to reduce network parameters and decompose the the 
network into subnetworks, if possible. If reliability is not completely computed via 
reliability-preserving reductions, the general SKT reliability expression is 
computed as described in chapter 4. If a specified CPU time limit is reached for 
any component, the computation is truncated and a lower bound on reliability is 
produced. 

There is a one main algorithm and numerous subroutine algorithms in the 
program. We discuss the main program and several important subroutina 
algorithms for understanding. The structure of the main program is based on a 
step by step approach which avoids redundant work. Comments will be given in 
curly brackets. 
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Main 



Input: A directed acyclic network D with vertex set V, V >2, edge set E, E =j 2, 
and set K<zV, \K\>2. Edge reliability Pi for each edge e t in E. 

Output: Rsk(D), exact reliability, if CPU time less than specified amount, otherwise 
a lower bound. 

(1) {Initialize} M = 1.0, bound=. false. 

(2) Read all edge data and create netviork data structure 

/ 

(3) {Change neck vertices to K-vertices} K = Neck ( D,K) 

(4) {Perform reductions for whole network} Reduct ( D.K.M ) • 

> 

(5) If remaining vertices | V\ < 1, go to end. 

(6) \D x ,Dz A} = Decomp(Z7) 

(7) For each subnetwork A 

(a) {Initialize} lowbnd=. false., set time = 0. 

(b) {Move edge tails to source} For each edge e =(i >i,Vj) with v^K 
{move Vi to s } Move(e ) 

(c) Construct the data structure for subnetwork A 

(d) {Reduce subnetwork} Reduct (A. K it M) 

(e) { Initialize } sumpro=0 

(f) { Find SAT reliability } sumpro = Rsa(A). 

(g) Count the remaining non- K-vertices for A. i.e., if vertex number 
> 0, count=count-l 




\ v U' v iz v ii> ~ CombO' ,i) 

2) { Find Hg, for £,VV-.;, . } M' = . . . ,v u ) 

3) sumpro =sumpro + M' • Rsa(ZJ— 

4) If elapsed CPU time & time limit then !owbnd=.true. and 
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bound=.true. 



5) If lowbnd=. false, and all combination have not been 
enumerated then go to (B.h. l) 

(i) M = M x sumpro 

(8) If bound=. false, print "Exact reliability is" M, otherwise print "Lower 
bound on reliability is" M 

End Main 

function Neck ( D,K) 

Input: D,K 

Output: New set of IC-vertices 

( This routine finds all neck vertices and changes them to K-vertices ) 

(1) Put all K into Q n 

(2) k = \K\ 

(3) While 0.^0 

(a) Remove v from Q n 

(b) Search all K-vertices in D-Vi using Depth-First Search 

(c) fc'= number of K-vertices reached 
(c) If k'<k then K=K[j\v J 

(3) Return (K). 

Eiid of Neck 



subroutine Reduct^.A",^) 

Input: D,K,M 

Output: Reduced network D with modified K and M 

( This routine performs all non-K-vertex reductions, series-parallel 
reductions and indegree-one K-vertex contractions ) 

(1) Construct Q d = $ veV—K \ outdeg(u)=0 J 
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(2) While Q d 

(a) For each vertex Vj incident to v it if outdeg(vj) = l and VjtV-K 
then add Vj to Q d . 

(b) Nonkre ( v t ) 

(3) Construct Q 3 = { v i e.V— K | indeg(^)= outdeg(i>i)=l } and 9c = \ V\tK \ 

indeg(^) = l I 

(4) For each if* 

(a) Check all outgoing edges, if e (v it Vj),e (vi,v } ) are found then 
Parall( e ) 

(b) If parallel reductions creates new series or contraction vertices, 
put these into 9, or 9c . respectively. 

(5) Yfhile 9s and 9c 

(a) If 9* *0 Series( Q,(i)) 

1) If series reductions create new parallel edges, do parallel 
reduction. Parall(e) 

2) If parallel reductions create new series or contraction 
vertices, put these into 9 S or Q c , respectively. 

(b) If 9 C *0 Contra( 9; (i)) 

1) If the contraction creates a new contraction vertex, put 
these into 9c 

2) If the contraction creates new parallel edges, do parallel 
reduction Parall^) 

3) If parallel reductions create new series or contraction 
vertices, put these into 9s or Q~, respectively 

(6) Return ( D,K,M ) 

End Redact 
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subroutine Decomp(Z?) 

Input: D 

Output: Subnetworks D\,D Z D m 

(1) While not all cut-vertices have been found 

(a) Using DFS, find pendant component Di with cut-vertex 

(b) If Di has no K-vertices except possibly then delete Di 

(2) Return (D h D z D m ) 

End Decomp 

After reading in the data and creating the network data structure, the 
algorithm begins by changing any neck vertices to K-vertices by a call to Neck at 
step (3). Next, a call to Reduct is made at step (4). Reduct first removes all 
outdegree zero non-K-vertices. Through the queue mechanism, it is ensured that 
all initial outdegree zero non-K-vertices are deleted along ■with any which are 
created as the reduction proceeds. A similar method is used for the other 
reductions schemes. 

The other reductions, indegree-one-K-vertex contractions and parallel and 
series reductions, may recursively create new reducible vertices or edges. The 
newly created reducible vertices are put into a queue instead of immediately being 
reduced to avoid redundant work. Parallel edges are reduced immediately, 
however. The network of Figure 5.1 is completely reduced by reductions alone and 
illustrates how a reduction can create other reducible edges or vertices. A call to 
Reduct would result in the set of reductions given below. This network would be 
completely reduced and its reliability completely calculated by these reductions. 
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(1) Reduce parallel edges 


e 3. e 4 to e 3 


(2) contract vertex 


v e 


( 3 ) series vertex 




( 4 ) contract vertex 


v 2 


( 5 ) reduce parallel edges 


2s, 25 to e 3 


(6) series vertex 


v 7 


(6) series vertex 


Vy 


( 7 ) reduce parallel edges 


e io, e ii to e 10 


(8) series vertex 


V 8 


( 9 ) reduce parallel edges 


2 10,2 12 to e 10 


(10) contract vertex 


^9 


(11) reduce parallel edges 


e 2 ,SQ to e 2 


(12) contract vertex 


v 3 


( 13 ) reduce parallel edges 


ej.Sy to 2 3 


( 14 ) contract vertex 


V 5 



After reductions are performed, a network will decompose into 
subcomponents if it contains any cut -vertices. This is done at step (6) of Main by a 
call to Decomp. The rest of the calculations Eire performed iteratively for each 
component A under step ( 7 ). First, all edges with tail-vertices in K are moved to 
the source. This may produce additional reductions, so another call to Reduct is 
made for the subcomponent. Then, the terms corresponding to Equation 4.2 (for 
that component) are produced until all such terms have been enumerated or the 
CPU time limit is reached. Note that this CPU limit is based on the CPU time used 
in calculating subcomponent reliability and is not based on total CPU time. 

If the CPU time limit is reached for any of the subcomponents, the Anal value 
obtained by the algorithm is a lower bound on reliability. Otherwise, the value 
obtained is exact reliability. Choosing a CPU time limit will depend on computer 
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facilities available and the network under analysis. Figure 6.5 shows the reliability 
increment as a function of CPU time. The lower bound obtained using a reasonable 
CPU time limit will not differ much from the exact reliability for this large scale 
network example 

B. DATA STRUCTURE 

There a many ways to represent a network on a digital computer. The most 
straightforward is to use an adjacency matrix: If vertices i and Vj are adjacent. 
M(i,j)=l, and M(i,j)=Q otherwise. Matrices of the above type are, in practice, 
extremely sparse and this representation is very inefficient in terms of space. It is 
often inefficient in terms of execution time, too, particularly when one is 
interested in operations that require retrieving information from all vertices 
adjacent to a particular vertex, e.g., Neck, Decomp, Reduct, and Rsa. In this 
situation, since most of the elements of M(i,j) are null, the algorithm will spend a 
great amount of time retrieving and comparing zero values. For a static network, 
an efficient representation is a packed matrix. In our algorithm, the network is 
dynamic, since it is reduced and decomposed many times during execution. 
Therefore, a multi-linked data structure (multiple linked lists with links to other 
data structures such as arrays) is used to represent the network. Every vertex has 
a Linked list of adjacent vertices, which, besides informing which vertices are 
adjacent to it, also tells whether or not the vertex belongs to set K. Two such lists 
are kept for each vertex, an "adjacent out” list and an "adjacent into" list. 0! 
course, each adjacency entry corresponds to an edge. So, in addition, there are 
pointers indicating the addresses where information about edges is kept. 

Figure 5.3 illustrates this data structure. The vertex only can send 
communication to vertices v 2 , and v 2 , and vertex v 2 can receive the 
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communication from vertex iq and send the communication to vertices v 3 and v 4 , 
and so on This structure is a very efficient data structure for dynamic networks. 

C. ALGORITHMIC COMPLEXITY 

It is important to determine the complexity of the various parts cf the 
algorithm so that overall complexity can be understood. In this section, we analyze 
the worst-case complexity of the subroutines, as written, and thus, find the 
complexity of the overall algorithm. Not all of the subroutines have been written 
as efficiently as possible because cf the difficulty in programming such routines 
and because of limited time. Furthermore, for problems of moderate size, it is 
unlikely that much efficiency would be gained using more sophisticated routines, 
since, in practice, most of the routines written seem to operate mere efficiently 
than the worst-case analysis indicates. 

There are three different reduction schemes included in the algorithm: non- 
K-vertex deletion, indegree-one K-vertex contraction and series and parallel 
reductions. However, the last three of these interact directly. 

The non-K-vertex deletion steps (1) and (2) in REDUCT, requires 0(|E’|) 
operations in the worst case. For example, consider a complete acyclic graph 
which has only one K-vertex, the source. Here, one vertex at a time would be 
deleted from the network, but, in the process, every edge would be examined 
exactly once until a single isolated vertex remained. The initial queue -building is cf 
no consequence since it is 0( | V\). 

Steps (3) through (5) constitute the other three reductions. Queue-building 
again may be disregarded since it is at an 0(1 V\) operation. Initial parallel 
reductions may require 0(|£’|) operations using a bucket sort technique to check 
for edges with common end vertices. The central step of the algorithm, step (5), 
requires 0(|7| 2 ) time. This is true since (1) at most | VJ series reductions or 



30 



contractions can ever occur (2) each such reduction requires one 0( K|) check for 
a newly created parallel edge, and (3) once identified, actually carrying out the 
parallel reduction requires only constant time. Overall then. REDUCT requires 
0( | K| a ) time. 

There are two different decomposition schemes used, cut-vertex 
decomposition and neck-vertex decomposition. The complexity of the cut-vertex 
algorithm Is 0(|£’|) since it is based on the 0(|£’|) depth-first search. The 
complexity of the neck-vertex algorithm will be 0( | V—K | |fT]) since finding neck- 
vertices is based on repeated 0( IfTl) searches on D-v for each non-K-vertex v . In 
the worst case then, this decomposition requires 0( | V] | £” | ) time. However, if we 
fix the number of non-K-vertices as described in chapter 4, then the complexity of 
neck-vertex algorithm will effectively be 0( |£’|). 

It follows from the above discussion that the total complexity of the reduction 
and decomposition algorithms is 0(! T /||iT|) normally, and OdFl 2 ) when non-K- 
vertices are fixed. 

After reductions, computing SKT reliability requires 0( \ E \2< V ~ K \) time. This is 
obviously true since combinations of non-K-vertices must be produced and 

an 0(| ill) SAT computation made for each combination. T iVhen the number of non- 
K-vertices is fixed, the computation becomes 0(|2T|). The entire algorithm 
therefore requires 0{\E\2^ V ~ K ') time in general but only 0(| V\ 2 ) time with non-K- 
vertices fixed. 
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VI. COMPUTATIONAL 
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In order to test the efficiency of the SKT reliability algorithm, we used two 
kinds of networks, "complete*' acyclic networks, and acyclic "street” networks. A 
complete acyclic network is a graph D(V,E ) such that V=\v l ,v z , • • • ,v n \ and 
e =(vi,Vj)e.E if and only if i<j. Any such network with n vertices has exactly 

71 edges. A network is called an mxn "street" network if (l) its vertices can 

& 

be laid out in an m. rows by n columns rectanguiar grid, (2) for all vertices except 
those in the last column there is an edge incident out each vertex directed to the 
right, and (3) for all vertices except those in the first row, there is an edge incident 
out of each vertex directed upward. When n=m. case the street network is called 
"square" (See Figure 6.1). 

In the complete network w r e assigned to be the source and fixed v n as a K- 
vertex. This insures that all edges and vertices cure relevant. Other X-vertices were 
chosen randomly as required. In the street network, all corner vertices were fixed 
as K-vertices, the lower left-hand corner vertex assigned as the source and all 
other K-vertices chosen randomly. 

A. S OUKCE-TO- ALL-TERMINAL RELIABILITY COMPUTATION 

As described in chapter 4. the complexity of computing SAT reliability is 
0(|£’|). Since repeated computation of SAT reliability is the key to computing 
general SKT reliability, it is interesting to see how much CPU time is required for a 
single SAT computation. Figure 6.2 shows the CPU time used as a function of the 
number of vertices in the two types of test networks. | V\ is increased from 15 to 
80 for the complete network and from 9 to 144 for the square street network. Of 
course, as a function of \V\, |£*| is quadratic for the complete network ranging 
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VI. COMPUTATIONAL EXPERIENCE 



In order to test the efficiency of the SKT reliability algorithm, we used two 
kinds of networks, "complete" acyclic networks, and acyclic "street” networks. A 
complete acyclic network is a graph D{V,E) such that V=\v u v 2 , ■ ■ ■ /u n } and 
e =(vi ,Vj)zE if and only if i<j. Any such network with n vertices has exactly 

n ( 7 ^~ 1 ) - edges. A network is called an mx.n "street" network if (l) its vertices can 

be laid out in an m rows by n columns rectangular grid, (2) for all vertices except 
those in the last column there is an edge incident out each vertex directed to the 
right, and (3) for all vertices except those in the first row, there is an edge incident 
out of each vertex directed upward. When 71=771 case the street network is called 
"square" (See Figure 6.1). 

In the complete network we assigned v x to be the source and fixed v n as a K- 
vertex. This insures that all edges and vertices are relevant. Other X-vertices were 
chosen randomly as required. In the street network, all corner vertices were fixed 
as K-vertices, the lower left-hand corner vertex assigned as the source and all 
other K-vertices chosen randomly. 

A. SOURCE-TO- ALL-TERMINAL RELIABILITY COMPUTATION 

As described in chapter 4, the complexity of computing SA*T reliability is 
0(|£'|). Since repeated computation of SAT reliability is the key to computing 
general SKT reliability, it is interesting to see how much CPU time is required for a 
single SAT computation. Figure 6.2 shows the CPU time used as a function of the 
number of vertices in the two types of test networks. | V\ is increased from 15 to 
60 for the complete network and from 9 to 144 for the square street network. Of 
course, as a function of \V\, | S’ ! is quadratic for the complete network ranging 
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from 105 to 3160. j^l increases linearly as a function of \V\ for the street 
network r anging from 12 to 284. From the figure, it is clear that SAT reliability can 
be computed very rapidly. 

B. SOURCE -TO-K-TERM INAL COMPUTATION 

Here we investigate the work required by the algorithm for computing general 
SKT reliability and the efficiency gained by using reductions and decomposition. 
Initial reductions and decomposition require 0(| V\ j^l) time. After reductions and 
decomposition, the complexity of computing SKT reliability is 

I y I 

0(max$ '2' 1 ^ 1 1 i }) where C is the set of separable components of D and where 

ieC 

Vi~Ki is the set of non-K-vertices in component i. Consequently, we expect an 
exponential growth in the worst case for SKT computations. Figure 6.3 shows the 
growth in CPU time as the number of non-K-vertices is increased, with and without 
reductions and decomposition 

For the complete network, the number of non-K-vertices ranged from 3 to 14 
out of a total 20 vertices. The reductions and decomposition do not reduce 
computation much here since the complete network is so densely connected. On 
the other hand, computation times are significantly reduced when reductions and 
decomposition are used in the street network, primarily because many non-K- 
vertices are changed to K-vertices. Notice that the computation time required by 
the street network without reductions and decomposition is almost the same as for 
the complete network. This shows how complexity is dependent mostly on the 
number of non-K-vertices. 

C. COMPUTATION WITH NON-K-VERTICES FIXED 

If the number of non-K-vertices is fixed, i.e. \ V-K\=c for some constant e, 
then the complexity of computing SKT reliability is 0(| F| 2 )+0(c |£’|) - 0(2° ,£’|) 
which is OdVl 3 ). Thus, the complexity is polynomial instead of exponential. Of 
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course, if reductions and decomposition are not used, the complexity becomes 
0( | £" | ) which is better than 0(|7| 2 ). However, in practice, the reductions and 
decomposition are very efficient and do not seem to increase computation time. 
Figure 6.4 shows the increase in computation for both test networks as j V\ is 
increased. 

D. LOWER BOUND COMPUTATION 

As described in chapters 4 and 5, a lower bound on SKT reliability is obtained 
when the normal SKT reliability computation is truncated. Two heuristic 
arguments were given in chapter 4 indicating that this lower bound should be fairly 
good. Figure 6.5 shows, for actual computations, how the lower bound approaches 
the exact network reliability very quickly. The test networks were assigned 
different values of p t to show that the accuracy of the bound is not much affected 
by Pi. It appears that a CPU time limit can be used successfully for truncating SKT 
computation to obtain a good lower bound when exact computation is not feasible. 
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VII. CONCLUSION 



This thesis has explored the source-to-K-terminal reliability analysis problem 
in acyclic directed networks. Vfe have shown that the polynomial-time reductions 
and decompositions reduce an exponentially complex problem to a polynomially 
complex problem under certain conditions. Even when these conditions do not 
hold, surprisingly large networks can be analyzed. The theoretical and 
computational results obtained here can be used by researchers in analyzing the 
reliability of directed acyclic networks and, in some cases, for the partial 
reliability analysis of general directed networks. 

“We assumed that all vertices work perfectly and that edge reliabilities are 
constant while the network system is alive. These assumptions may not be 
acceptable in a real-world situation. It may be necessary to treat vertex 
reliabilities explicitly and. in some way, handle component reliabilities changing as 
a function of time. This will complicate reliability analysis but does not make 
analysis infeasible. 

The extension to unreliable vertices is not hard. All K-vertices must work if 
the system is to work and thus, for unreliable K-vertices, we obtain network 
reliability as usual and then multiply by the probability that all K-vertices function. 
Unreliable non-K-vertices may be handled in several ways. One of the simplest 
ways is to split each unreliable non-K-vertex into two perfectly reliable vertices 
with an unreliable edge between the two vertices with a failure probability equal to 
the failure probability of the original vertex. All edges incident into the original 
vertex are made incident into the tail vertex of the new edge and all edges incident 
out of the original vertex are made incident out of the new head vertex. However. 
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it is also possible to redefine the reduction and decomposition schemes and 
redefine the general reliability expression to directly handle unreliable non-K- 
vertices. For example, the series reduction can include the vertex reliability so 
that p c = Also, the term in the general reliability expression associated 

with "no edges in Ey..* work" can be replaced with. a term which expresses "no 
edges into -v* work or v. : fails, and no edges into Vj work or v, faiis, and ... except 
edges going between pairs of vertices in [Vi.Vj, • • • ,v k l. This is not hard to 
compute. 

Reliability changing as a function of time can be handled in one of two ways. 
The easiest way is to just input component reliabilities at different times and 
compute network reliability at these times. Interpolation can be used to 
approximate values between the times for which reliability is actually evaluated. 
Another possibility is to produce a symbolic reliability expression that can be 
repeatedly evaluated by a computer at different times, thus avoiding repeated 
reductions, decompositions, etc. This might require large amounts of storage but, 
if reliability can be calculated at all in a reasonable amount of time, it should be 
feasible. Once the price is paid, computing reliability for any given point in time 
should be very efficient. 

As modern computer, power and other network systems become larger and 
more complicated, reliability becomes a more important part of design criteria. 
Unfortunately, it is obvious that as networks become more complicated they also 
become more difficult to analyze. Further research should concentrate on ways to 
streamline computations and to find accurate but easily computable bounds. The 
method described in chapter 4 for computing SKT reliability could be improved if 
only non-zero terms were generated; An enumeration scheme similar to that of 
Satyanaryana and Prabhakar [1978] should help in this respect. The lower bound 

t 

described in this thesis is a good start on finding accurate and easily 'computable 
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bounds, but every lower bound should be accompanied with an upper bound so that 
accuracy can be checked. Finding a good upper bound should be a priority for 
future research. 
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Neck Vertex to K-vertex Reduction 
< Figure 3.1> 





Decomposition 
<Figure 3.2> 
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Moving Edges to Source 
<Figure 3.3> 




All Terminal Reliability R s v{D\) = R s y(Dn) 
<Figure 4.1> 
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Completely Reducible Network by Reduction Schemes 

< Figure 5. 1 > 




Edge Vertex 




AdJ . In V AdJ , out 



Multilist Data Structure 
<?ig ore 5.2 > 
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4x4 Square "Street” Network 
<Figure 6.1> 
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Source-To-All Terminal Network 
< Figure 6.2 > 
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ELAPSED CPU TIC WITH ANO WITHOUT REDUCDON, OCCOMPQSffiON 




LEGEND 
o a SO. ST. WITH 
o s COMPLETE WTTH 
a * SQ.ST. WITHOUT 
♦» COMPLETE WITHOUT 



Soucs-To-if Terminal Network 



< Figure 6.3 > 
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COMPUTER CPU TIC WITH FIXED NON-K-VERT1CE5 10 




LEGEND 

o ■ COMPLETE fCTWORK 
o * STREET SO. NETWORK 



Score e-To-/T Terminal with Fixed Non- A"-' Vertices 
< Figure 6.4 > 
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RQLAMJTY WCREASE BY TIME IN TOTAL 20 VERTICES 
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Network Reliability Increment 
< Figure 6.5 > 
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